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ABSTRACT 

Simultaneous  measurements  of  wave  height  and  two  orthog- 
onal water  particle  velocities  were  made  at  six  elevations 
in  nineteen  meters  of  water  using  a  penetrating  wave  staff 
and  an  electromagnetic  flowmeter.   Moderate  swell  and  low 
wind  conditions  prevailed  during  the  experiment.   The  mea- 
sured wave-induced  velocities  were  two  to  four  percent 
greater  than  those  predicted  using  linear  wave  theory. 
Coherence  of  the  wave  height  and  wave-induced  velocities  in 
the  significant  energy-density  range  was  computed  to  be  over 
0.85,  indicating  that  the  motion  was  almost  totally  wave 
induced.   At  higher  frequencies  it  was  apparent  that  the 
motion  was  primarily  turbulence.   Phase  spectra  computed  for 
the  measured  wave  heights  and  orbital  velocities  compared 
very  well  with  linear  theory.   Measured  frequency  distribu- 
tions were  compared  to  both  Gaussian  and  Gram-Charlier 
distributions  using  the  chi-squared  goodness-of-fit  test. 
Qualitatively,  the  Gram-Charlier  distribution  gave  the 
better  fit  to  the  data. 
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I..   INTRODUCTION 

A.   BACKGROUND 

Measuring  instantaneous  water  particle  velocities  in 
the  presence  of  wave-induced  motion  has  been  difficult  due, 
primarily,  to  a  lack  of  instrumentation.   This  is  particu- 
larly true  for  measurements  made  in  the  field.   Until  re- 
cently, there  has  been  no  adequate  instrument  available  to 
measure  rapidly  fluctuating  flows  found  beneath  wind  waves 
and  swell.   As  well  as  the  usual  requirements  of  reliability 
and  accuracy,  such  an  instrument  must  also  be  rugged,  have  a 
low  response  time,  and  have  good  directional  measuring 
characteristics . 

Earlier  attempts  at  correlating  water  height  and  particle 
motion  have  generally  been  accomplished  in  shallow  water  next 
to  the  shore,  off  piers,  or  in  estuaries.   Historically, 
Inman  and  Nasu  (1956)  made  one  of  the  earliest  attempts  at 
field  measurements  of  wave-induced  water  particle  motion, 
using  a  force  meter  in  shallow  water  and  compared  their 
results  to  solitary  wave  theory.   Miller  and  Zeigler  (1964) 
used  both  an  acoustic  meter  (based  on  the  doppler  shift 
principle)  and  an  electromagnetic  flowmeter  to  measure  water 
particle  motion  in  the  surf  zone.   They  compared  their 
measurements  to  higher  order  wave  theory  and  obtained 
qualitative  agreement. 

Nagata  (1963,  1964a,  1964b)  used  an  electromagnetic 
flowmeter  in  a  series  of  experiments  to  measure  water 


10 


particle  velocities  in  the  surf  zone  and  to  determine  wave 
directional  spectra.   Bowden  and  White  (1966)  and  Simpson 
(1968)  performed  field  measurements  using  an  electromagnetic 
flowmeter  developed  at  the  National  Institute  of  Oceanography 
Like  Nagata,  Bowden  and  White  (1966)  used  water  particle 
velocities  to  determine  the  directional  spectrum  of  waves. 

Simpson's  experiment  was  performed  at  an  open  coastal 
site  in  water  approximately  6  meters  deep.   The  waves, 
primarily  wind  generated  sea  waves,  were  developed  over  a 
limiting  fetch  of  210  km.   Simpson  employed  spectral  methods 
in  the  analysis  of  the  waves  and  of  the  wave-induced  water 
particle  velocities.   His  results  were  consistent  with 
linear  wave  theory. 

Seitz  (1971)  carried  out  a  similar  experiment  in  that  he 
measured  sea  surface  displacement  and  water  velocity  compo- 
nents in  an  estuary  at  a  depth  of  one  meter.   He  compared 
the  measured  velocity  spectra  to  a  derived  velocity  spectra 
obtained  from  direct  measurement  of  the  surface  wave  field 
using  linear  theory.   Seitz  had  to  contend  with  a  shift  in 
the  observed  frequencies,  the  magnitude  of  which  was  a 
function  of  the  tide  induced  velocity.   Since  at  the  fre- 
quency of  peak  amplitude  the  error  introduced  was  approxi- 
mately 6.7%,  he  did  not  correct  the  data.   The  total  water 
particle  motion  at  depth  can  be  described  as  the  sum  of  the 
turbulent  velocities  and  the  surface  wave  induced  motion. 
Seitz  attempted  to  separate  the  two  based  on  the  slope  of 
the  spectrum  in  the  inertial  sub-range. 


11 


The  lack  of  instrumentation  has  not  been  the  only 
hindrance  to  v:ave  height  and  velocity  studies.   Most  research 
has  been  conducted  at  shallow  water  inshore  sites  because  of 
the  general  unavailability  of  a  deep  water,  exposed  instru- 
ment platform.   Shonting  (1966)  used  the  Buzzard's  Bay 
Entrance  Light  Station  on  the  South  Coast  of  Massachusetts 
to  obtain  a  twenty  meter  water  depth  at  an  open  ocean  site. 
He  used  ducted  current  meters  to  measure  the  orthogonal 
components  of  the  water  velocity  under  waves.   His  experiment 
was  designed  primarily  to  evaluate  these  meters.   An  inherent 
difficulty  with  any  type  of  ducted  meter  is  their  poor 
directional  measuring  characteristics  when  the  angle  between 
meter  axis  and  the  flow  is  large. 

Bordy  (1971)  made  measurements  in  20  meters  of  water 
from  the  Naval  Undersea  Research  and  Development  Center  (NUC) 
Oceanographic  Research  Tower  located  approximately  one  mile 
off  Mission  Beach,  San  Diego.   He  compared  wave  height  and 
water  particle  velocity  using  spectral  analysis,  and  found 
a  general  agreement  with  linear  wave  theory  except  for 
measurements  of  phase  spectra  at  the  deepest  measurements  of 
ten  meters.   The  phase  shift  encountered  was  unexplained  in 
the  phase  angle  error  analysis.   During  the  experiment,  the 
NUC  tower's  permanently  installed  wave  staff  system  was  used 
for  water  level  sensing.   The  flowmeter  was  removed  horizon- 
tally from  the  wave  staff  by  approximately  five  feet.   It  is 
felt  that  due  to  this  displacement,  an  error  was  introduced. 
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During  8  June  1972  records  of  orthogonal  water  particle 
velocity  and  instantaneous  water  elevation  were  again 
measured   at  the  NUC  tower.   The  data  were  taken  in  conjunc- 
tion with  an  experiment  involving  the  investigation  of  wave 
induced  fluctuations  in  acoustic  phase  shifts,  temperature, 
and  salinity. 

The  results  of  the  wave  and  water  particle  motion 
measurements  are  described  here.   An  electromagnetic  flow- 
meter was  used  to  gather  the  wave  particle  velocity  data 
at  several  depths.   These  measured  values  of  wave-induced 
motion  were  compared  with  values  predicted  using  linear 
theory  and  the  wave  height  record. 

B.   RESEARCH  OBJECTIVE 

The  purpose  of  this  investigation  was  four  fold: 

1.  To  measure  simultaneously  wave  height  and  two 
orthogonal  wave-induced  water  particle  velocities. 

2.  To  perform  a  spectral  analysis  on  the  data  to  obtain 
its  statistical  properties. 

3.  To  spectrally  compare  the  data  to  first  order  wave 
theory. 

4.  To  determine  the  statistical  geometry  of  both  the 
sea  surface  and  the  measured  velocities. 
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II.   EXPERIMENTAL  SITE  AND  INSTRUMENTATION 

A.  EXPERIMENTAL  SITE 

The  experiment  was  conducted  at  the  Naval  Undersea 
Research  and  Development  Center  Oceanographic  Research  Tower 
at  San  Diego,  California,  8  and  9  June  1972.   The  NUC  tower 
is  located  approximately  one  mile  off  Mission  Beach,  San 
Diego  in  a  mean  water  depth  of  19  meters.   It  is  constructed 
of  steel  and  concrete  and  is  permanently  imbedded  in  the  sea 
floor,  Figure  (1).   The  forty  foot  square  tower  has  three 
levels.   The  first  is  used  primarily  for  gaining  access, 
while  the  second  houses  machinery  and  provides  work  space 
for  placing  instruments  in  the  water.   The  third,  which  is 
enclosed,  is  used  for  the  housing  of  electrical  equipment 
and  also  provides  sleeping  space  for  six  men.   The  legs  of 
the  tower  slope  five  degrees.   A  pair  of  tracks  is  mounted 
on  the  north,  south,  and  west  sides  for  lowering  scientific 
instrument  packages  into  the  sea  to  any  depth  above  the 
ocean  floor.   The  western  side  of  the  tower  was  used  for  the 
experiment  because  of  its  exposure  to  the  predominant  swell. 

i 

B.  INSTRUMENTATION 

Water  particle  velocities  were  measured  with  an  Engineer- 
ing Physics  Company  water  current  meter  Model  EMCM-3B.   The 
sea  surface  elevation  was  measured  with  a  Baylor  Company 
wave  staff  system  Model  13528R.   The  penetrating  wave  mea- 
suring system  gives  a  direct  measure  of  the  sea  surface 
elevation. 
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Figure  1.   Naval  Undersea  Research  and  Development  Center 
Oceanographic  Research  Tower,  San  Diego,  California 
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The  model  13528R  Baylor  Company  V/ave  Staff  System  is  an 
instantaneous  water  level  measuring  device  consisting  of 
two  tensioned   1/2",  6x19  IWRC  304  SS  wire  ropes.   Together 
with  these  lengths  of  stainless  steel  wire  rope,  the  trans- 
ducer produces  an  electrically  linear  direct  current  output 
that  is  proportional  to  the  amount  of  v/ave  staff  above  a 
short  circuit  produced  by  the  water  surface.   System  accuracy 
is  stated  as  within  one  percent  of  actual  height  and  the  time 
constant  is  less  than  0.06  seconds. 

As  received  from  the  manufacturers,  the  working  length 
of  the  staff  was  50  feet;  however,  in  order  to  take  water 
particle  velocity  measurements  near  the  surface  directly 
under  the  staff,  it  was  shortened  to  20  feet.   After 
shortening,  the  staff  was  recalibrated.   The  resulting 
calibration  curve  is  shown  in  Figure  (2).   The  rails  on 
which  the  instrument  cart  rides  slopes  with  the  legs  of  the 
tower.   This  necessitated  the  rigging  of  an  aluminum  "I" 
beam  with  a  trolley  to  support  the  Baylor  wave  staff  system 
directly  over  the  electromagnetic  flowmeter.   The  wave  staff 
was  tensioned  by  allowing  the  weight  of  the  instrument  cart 
and  frame  to  be  exerted  on  the  lower  end  of  the  wave  staff 
cables.   The  estimated  tension  was  approximately  800  pounds. 
The  instrument  and  installation  are  shown  in  Figure  (3). 

The  EMCM-3B  electromagnetic  current  measuring  device 
measures  two  orthogonal  components  of  water  particle  velocity 
through  a  range  of  zero  to  five  M/SEC.   Other  characteristics 
of  the  flow  meter  are: 
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Figure  3.   Model  13528R  Baylor  Company  Wave  Staff  Syste 
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Output  Voltage: 
Electrical  Time  Constant 

Error  Band: 


Transducer  Size: 


±5  volts  full  scale 

0.2  seconds  or  longer, 
adjustable  by  user 

the  largest  of 

1)  ±1  percent  of  full  scale 

2)  ±5  MM/SEC 

3)  5/(T^)  MM/SEC 

where  T  is  the  time  constant 

2  3/4"  O.D.  x  36"  L 


The  operation  of  the  unit  is  based  on  Faraday's  principle  of 
magnetic  induction.   A  current  is  induced  when  a  conducting 
fluid  flows  around  the  probe  through  an  alternating  magnetic 
field  generated  by  a  toroidal  magnet  inside  the  instrument. 
The  equation  describing  the  production  of  current  is: 

E  =  u  X  S, 

where  E  is  the  induced  voltage,  u  the  velocity  of  the 
conducting  medium,  and  B  the  magnetic  field.   The  induced 
electromotive  force  is  linearly  proportional  to  the  amount 
of  fluid  flowing  through  the  magnetic  field.   The  magnetic 
field,  B,  is  most  intense  in  the  near  vicinity  of  the  probe 
and  decreases  according  to  the  inverse  square  of  the  distance 
from  the  instrument.   Accordingly,  the  flow  velocity  near 
the  probe  is  weighted  most  heavily  in  its  contribution  to  the 
measured  velocity.   Qualitatively,  only  that  water  which  is 
within  a  distance  of  two  to  three  probe  radii  contributes 
significantly  to  the  flow  induced  voltage. 
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The  flow  meter  was  originally  calibrated  under  steady 

flow  conditions  by  the  manufacturer.   Different  response 

characteristics  are  to  be  expected  under  unsteady  flow 

conditions;  consequently,  the  probe  was  recalibrated  by 

oscillating  it  in  a  water  tank.   The  equipment  schematic  is 

shown  in  Figure  (4).   The  carriage  on  which  the  flow  meter 

probe  was  mounted  travels  back  and  forth  on  rails.   It  is 

driven  by  a  variable  speed  electrical  motor.   The  peak 

carriage  velocity  was  calculated  from  the  tangential 

velocity  of  the  motor  arm.   The  ratio  of  carriage  velocity, 

V,  ,  to  the  velocity  measured  by  the  instrument,  V  ,  was 
l>  in 

found  for  different  angular  velocities.   This  was  done  for 
each  pair  of  electrodes  by  orienting  them  parallel  to  the 
flow  in  turn. 

The  instrument's  response  was  very  linear  with  respect 
to  frequency  up  to  at  least  0.5  Hz.   Beyond  this  frequency 
the  cart  upon  which  the  probe  was  mounted  became  unstable. 
Thus  0.5  Hz  was  the  maximum  frequency  investigated  during 
calibration.   This  limitation  was  not  deemed  a  debilitating 
factor  since  0.5  Hz  was  above  the  upper  frequency  limit  of 
significant  wave  energy.   The  results  of  the  calibration  are 
shown  in  Figure  (5)  and  the  instrument  and  installation  in 
Figure  (6). 
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Figure  5.   Calibration  Curve  for  the  Engineering  Physics 
Company  Model  EMCM-3B  Flowmeter 
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Figure  6.   Engineering  Physics  Company  Model  EMCM-3B 
Electromagnetic  Flowmeter  Installed  in  Position  (arrows) 
on  the  instrument  cart 
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III.   EXPERIMENTAL  PROCEDURE 

Sea  surface  elevation  and  two  orthogonal  components  of 
water  particle  velocity  were  measured  simultaneously.   The 
cart  on  which  the  flowmeter  was  attached  was  raised  and 
lowered  so  that  measurements  could  be  made  at  various  depths. 

Data  from  the  flowmeter  and  wave  gauge  were  recorded  on 
a  Sangarno  Model  3500  fourteen  channel  FM  tape  recorder 
during  the  ten  data  collection  periods.   The  data  chosen  for 
analysis  was  obtained  during  the  first  six  recording  periods 
on  8  June  1972.   Each  processed  run  was  about  20  minutes 
duration.   After  initial  processing,  the  other  four  data 
collection  periods  were  deemed  too  short  to  be  statistically 
representative . 

Summarized  in  Table  I  below  are  the  velocity  components, 
water  depth,  depth  of  the  instrument,  and  record  length  for 
each  of  the  processed  runs.   Instrument  depth  refers  to  the 
depth  of  the  flowmeter  below  mean  sea  level.   The  u-component 
of  velocity  was  oriented  north-south  approximately  parallel 
to  the  shore,  while  the  v-component  was  measured  east-west. 
Correct  flowmeter  orientation  was  obtained  using  a  plumb-bob 
and  level.   The  wave  gauge  was  always  maintained  directly 
over  the  flowmeter.   As  the  cart  was  lowered  the  top  of  the 
wave  staff  was  pushed  further  away  from  the  tower  to  compen- 
sate for  the  five  degree  tilt  of  the  tracks  to  which  the 
cart  was  fastened. 
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The  wave  staff's  perpendicularity  to  the  horizontal  was  also 
ascertained  by  the  use  of  a  level.  The  wave  staff  installa- 
tion is  pictured  in  Figures  (7)  and  (8). 


TABI 

,E  I 

RUN 

DEPTH 
(meters ) 

DEPTH  OF 
INSTRUMENT 
(meters ) 

COMPONENT 

OF 
VELOCITY 

LENGTH  OF 
RECORDS 
(minutes ) 

1 

19.15 

7.04 

u 

&  V 

24.96 

2 

19.22 

14.62 

u 

&  w 

22.61 

3 

19.25 

12.64 

u 

&  w 

21.33 

4 

19.28 

IO.69 

u 

&  w 

24.85 

5 

19.32 

8.66 

u 

&  w 

18.88 

6 

19-36 

5.70 

u 

&  w 

19.52 

The  sea  surface  was  in  general  disturbed  by  swell  corning 
from  the  west-southwest.   The  sky  was  overcast  with  a  low 
stratus  layer  and  little  wind,  less  than  10  knots  generally. 
Sea  surface  temperature  was  obtained  using  a  bucket  thermom- 
eter.  Table  II  summarizes  the  environmental  data. 

TABLE  II 

TIME  OF    AIR     SEA   \  WIND      CLOUD   SWELL 

RUN   START  OF   TEMP    SURF         DIR      VEL   COVER   DIR 

RUN  (PPT)  (°C)    TEMP  (°C)    (°T)   (knots)  (%) 


1 

1405 

19.7 

19.20 

240 

8 

100 

wsw 

2 

1600 

19.6 

18.95 

250 

10 

100 

wsw 

3 

1654 

19.2 

18.20 

240 

9 

100 

wsw 

4 

1800 

18.7 

19.10 

225 

6 

100 

wsw 

5 

1830 

18.2 

18.99 

225 

7 

100 

wsw 

6 

1922 

17.9 

18.77 

235 

4 

100 

wsw 
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Figure  7.   Model  13528R  Baylor  Wave  Staff  Installation 
Showing  "I"  Beam  and  Trolley  for  the  Top  of  the  System, 


m 


■ 


Figure  8.   Pulley  Attachment  Point  for  the  Bottom  of  the 
Baylor  Wave  Gauge  (arrow). 
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Table  III  summarizes  the  temperature  profile  data  taken 
at  the  beginning  of  each  data  collection  period.   This  data, 
which  shows  a  remarkable  variation  in  layer  depth  and 
intensity  of  temperature  gradient,  was  obtained  using  a 
mechanical  BT.   The  NUC  tower  was  originally  constructed  to 
study  internal  waves  which  are  prevalent  in  the  area.   The 
variability  in  temperature  profiles  is  attributed  to  internal 
waves.   A  typical  temperature  profile  is  shown  in  Figure  (9). 

TABLE  III 


RUN 

BT 
SURF 
TEMP ( °  C ) 

1 

22.8 

2 

22.5 

3 

23.3 

4 

21.1 

5 

22.2 

6 

19.9 

TOP  OF  BOTTOM  OF  TOP  OF  BOTTOM 

1st  LAYER  1st  LAYER  2nd  LAYER  TEMP 

(°C/m)  ("C/m) (°C/m) (°C) 

22.2/7.31  16.6/12.19  xxxxxxxxxx  15-55 

22.2/3.04  19.9/12.80  17.2/15.85  16.66 

23.3/6.08  21.6/9.1^  20.5/16.46  18.33 

21.1/4.87  19.4/6.71  18.3/17.68  16.66 

21.6/5.49  19.9/10.67  17.7/17.06  18.33 

19.9/5.79  16.1/15.24  xxxxxxxxxx  16.11 


The  analog  FM  recordings  were  digitized  using  a  hybrid 
system  consisting  of  a  Scientific  Data  Systems  (SDS)  Model 
5000  analog  computer  and  a  Xerox  Data  Systems  Model  9300 
digital  computer.   For  comparison  purposes  the  data  were 
transcribed  onto  a  rectilinear  eight  channel  Clevite  strip 
chart  recorder  after  a  tenfold  amplification  through  the 
(SDS)  5000  analog  computer.   An  example  of  the  transcribed 
data,  after  digitization,  is  shown  in  Figure  (10).   After 
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Figure  9.   Typical  Temperature-Depth  Profile 
Run  Number  4.   8  June  1972. 
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digitizing  onto  the  seven  track  tape  used  by  the  9300  system, 
the  data  were  transcribed  onto  a  nine  track  tape  using  the 
Naval  Postgraduate  School's  IBM  Model  360  digital  computer. 
This  computer  was  also  used  for  the  statistical  analysis. 
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IV.   THEORY 

A.   LINEAR  WAVE  THEORY 

The  elevation  of  the  surface  n(t)  can  be  described  as 
the  superposition  of  an  infinite  number  of  sinusoids  of  the 
form: 

CO  oo 

n(t)  =  I     a  cos  (£  -x-<j-  t+e  )   =  I     n        (4.1) 

n=l  n       n    n   n    n=l  n 

-»■ 

where  x  is  the  horizontal  Cartesian  coordinate,  t  is  the 

time,  k  is  a  horizontal  vector  wave  number,  a      the  frequency 

->- 
related  in  linear  theory  to  k  by: 

c  2  =  gk   tanh  k  h  (4.2) 

n    &  n       n 

where  g  is  the  acceleration  of  gravity,  a   is  the  amplitude 
and  e   the  phase  angle.   n(t)  represents  the  sum  total  of 
all  component  wavelets.   Summing  in  the  manner  of  Equation 
(4.1)  implies  a  linear  system  and  restricts  the  analysis  to 
the  use  of  linear  or  Airy  wave  theory  to  describe  the  wave 
induced  motion. 

Linear  wave  theory  can  be  used  to  predict  the  wave 
induced  water  particle  velocities.   The  equations  for  the 
horizontal  and  vertical  velocities  respectively  are: 


«  a^gk^  cosh  k^(h+z) 

cos(kn-x-ant+en) 

(4.3) 


u(t)  =  I     -?—^ 2 _  cos(kn.x-ant+en) 

n=l    a   sinh  k  h 
n       n 


00  gk  cosh  k  (h+z) 

r  n      n 

^•,  ~~o      sinh  k 

n=l  n       n 


r   "  n n_ 

A-,        a.  sinh  k  h     nn 
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and 


°°   angkn  sinh  kn^h+z)      ->   + 

w(t)  =   J   ; — r— i — i sin(k  'x-o   t+e  ) 

S     a   smb  k  h  n    n   n 

n=l     n       n 


00   gk   sinh  k  (h+z) 


(4.4) 


£     a 


n 


_-,    a   sinh  k  h 
n=l    n       n 


n 


-nno 


e=90 


where  h  is  the  total  v/ater  depth  and  z  is  the  depth  of 
interest  measured  positively  upward  from  the  still  water 
level. 

The  terms 


gk   cosh  k  (h+z) 
&  n n     ' 

a   sinh  k  h 
n       n 


(4.5) 


and 


gk   sinh  k  (h+z) 
&  n       n 


a   sinh  k  h 
n       n 


(4.6) 


represent  the  transfer  functions  relating  the  horizontal  and 
vertical  velocities  to  n(t). 

B.   LIMITATIONS  OF  FIRST  ORDER  WAVE  THEORY 

In  the  formulation  of  linear  wave  theory,  the  boundary 
conditions  are  linearized  in  order  to  obtain  an  analytical 
solution.   In  the  linearization,  it  is  assumed  that  the 
amplitude,  a,  is  small  compared  to  the  wave  length,  L,  that 
is,  a/L  <<  1. 

Higher  order  solutions  to  the  boundary  value  problem, 
generally  obtained  by  perturbation  analysis,  give  a  better 
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representation  of  a  constant  profile  wave.   However,  the 
nonlinearities  introduced  in  the  solution  precludes  their 
use  where  the  principle  of  superposition  is  invoked.   This 
in  turn  eliminates  their  generalization  to  a  stochastic 
process . 

Dean  (1968)  analyzed  the  various  wave  theories  with 
respect  to  how  well  they  fitted  both  the  dynamic  and  kine- 
matic free  surface  boundary  conditions.   He  showed  that 
first  order  theory  provided  a  good  solution  to  the  wave 
equation  for  waves  of  permanent  form  having  small  steepness 
in  intermediate  water  depths  (the  assumption  of  linear 
theory  stated  above). 

The  conditions  prevailing  during  this  experiment  fell 
within  these  limitations.   The  rrns  wave  height  averaged  0.27 
meters  with  a  mean  period  of  16  seconds  corresponding  to  a 
deepwater  wavelength  of  209  meters  in  19  meters  of  water. 
Hence,  if  linear  theory  provides  a  reasonable  approximation 
for  the  average  wave  conditions,  the  superposition  of  a 
large  number  of  components,  each  described  by  linear  theory, 
should  provide  an  even  better  description. 
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V.   THEORETICAL  CONSIDERATIONS  OF  DATA  PROCESSING 

BY  STATISTICAL  ANALYSIS 

The  digitized  records  were  spectrally  analyzed  using 
the  methods  of  Blackman  and  Tukey  (1958)  and  utilizing  the 
IBM  Model  360/67  computer  at  the  Naval  Postgraduate  School. 

A.   COMPUTATION  OF  THE  ENERGY-DENSITY  SPECTRUM 

The  covariance  function  d>      (t)  of  time  series  random 

xx 

data  describes  the  general  dependence  of  the  values  of  the 
data  at  one  time  on  the  values  at  another  time  for  the  same 
data  set.   In  integral  form 


*xx(t)  =  £S  T     I       x(t)  x(t+T)  dt  (3-1} 


which  is  the  average  of  the  product  of  the  value  of  a  data 
point  x(t)  with  that  of  an  adjacent  point  at  time  x  later. 
The  resulting  average  product  will  approach  the  exact  auto- 
correlation function  as  T  approaches  infinity. 

When  using  a  covariance  function  generated  with  a  finite 
data  set,  it  is  necessary  to  apply  a  lag  window  to  <$>  r    (t). 
A  Parzen  window,  as  used  in  this  analysis,  has  the  desirable 
property  of  no  negative  side  lobes  thus  eliminating  any 
possible  numerical  instability  which  might  arise  due  to 
resonance  effects  of  the  side  lobes.   This  is  particularly 
critical  in  computing  cross-spectra.   The  Parzen  window  is 
defined  as: 


3^ 


P(t)  =  1  -  6(x/Tm)2  +  6(i/Tm)3      t  =  0,1,2,. ..,Tm/2 
=  (1  -  (t/T  ))3  t  =  T/2  +  1,...,T 

=  0  t  -  T 


m 


(5.2) 


where  T   Is  the  total  lag  time, 
in  ° 


The  energy-density  spectrum  as  computed  from  the  Fourier 
transform  of  the  covariance  function  with  the  Parzen  window 
applied  to  correct  for  finite  record  length  is 

CO 

*xx(f)  =  I      P(T)  ^xx(t)  e-i27TfT  dx   ,  (5.3) 

_oo 

where  f  denotes  a  particular  frequency. 

The  energy-density  spectrum  is  the  statistical  estimate  of 
the  energy-density  distribution  in  the  various  frequency 
bands  comprising  the  spectrum  and  the  area  under  an  energy 
density  spectrum  is  equal  to  twice  the  variance  of  the  random 
variable.   With  respect  to  the  wave  height  spectrum3the  area 
is  proportional  to  the  total  potential  energy-density  of  the 
sea  surface  elevation,  while  the  area  under  the  velocity 
spectrum  is  proportional  to  the  total  kinetic  energy-density. 

B.   COMPUTATION  OF  THE  CROSS-SPECTRAL  DENSITY 

The  cross-covariance  function  is  similar  to  the  covariance 
fucntion  with  the  exception  that  it  compares  two  time  records 
lagged  by  a  time  interval  x.   Thus,  in  integral  form 
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lim   1 


,xy(T)  s  ?IZ     T-  I     x^)  y^t+T)  dt  •  (5.4) 


If  x  and  y  are  statistically  Independent  and  the  mean  of 

both  x  and  y  is  zero,  §      (x)  =  0  for  all  f.   If  not,  <j>   (t) 

xy  xy 

is  equal  to  the  product  of  the  means.   The  lag  window 
considerations  for  a  finite  record  length  are  identical  to 
those  for  the  covariance  function.   The  cross-spectral 
density  is  therefore  the  Fourier  transform  of  the  Parzen 
window  lagged  cross-covariance  function 

00 

*xy(f)  =  /   P(t)  d>xy(T)  e-i2TrfT  dx  (5.5) 

The  cross-spectrum  can  be  defined  in  terms  of  its  real 
and  imaginary  parts  since,  unlike  the  energy-density  spectrum, 
it  is  not  an  even  function.   Thus 


$   (f)  =  C   (f)  -  iQ   (f)  (5.6) 

xy       xy        xy 


where   the    real   part   of   $      (f)    is   the    co-spectrum 


C      (f)    =    2    /      [4      (t)    +   ♦__(-t)]    cos(27TfT)    dx    ,  (5.7) 

xy  q  Ay  Ay 


and  the  imaginary  part  the  quadrature  spectrum 


Q      (f)    =   2   /      C4>      (t)    -   <t>      (-t)]    sin(2irfT)    dx    .  (5.8) 

xy  q  -*-y  &y 
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C.   COMPUTATION  OF  COHERENCY 

A  useful  quantity  arising  out  of  the  computation  of  the 
cross-  and  energy-density  spectra  is  the  coherency. 
Coherency  is  defined  as 

Y2xv(f)  «  & ,  (5.9) 

Xy       $   (f)  $   (f) 

xx'    yy 


where 

- 


0  <  Y2_,(f)  <  1  .  (5.10) 


2 
When  y  Y_(f)  =  1,  for  all  f,  the  two  records  x  and  y  are 

xy 
fully  coherent  and  statistical  dependence  can  be  assumed. 

On  the  other  hand,  if  y   yv(f)  =  0  for  all  f,  x  and  y  are 

statistically  independent.   Bendat  and  Piersol  (1966) 

attribute  coherency  values  less  than  unity  to  three  possible 

causes:   the  system  is  not  linear,  extraneous  noise  is 

present  in  the  measurements,  or  the  output  response  is  due 

to  more  than  one  input  function. 

D.   COMPUTATION  OF  PHASE  ANGLE 

The  cross-spectral  phase  angle  represents  the  average 
angular  difference  by  which  the  cross-correlated  components 
of  y(t)  lead  those  of  x(t)  in  each  spectral  frequency  band. 
The  phase  angle  can  be  computed  using  the  co-  and  quadrature- 
spectra 

e   (f)  =  arctan  [Qxy(f)  /  cxy(f)J   •  C5.ll) 
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E.   PROBABILITY  DENSITY  FUNCTION  OF  THE  SEA  SURFACE 
AND  INDUCED  WATER  PARTICLE  FLUCTUATIONS 

The  probability  density  function,  p(x),  of  random  data 

describes  the  probability  that  the  data  will  assume  a  value 

within  some  defined  range  at  any  instant  of  time. 

Mathematically 

,  v  _  lim   Prob[x<x(t  )<(x+dx)  ]  _  lim   JL_  rlirn  ^x_-,    , ^   _p* 
PKX)        dx+0  dx  dx-^0  dx  LT+«>   TJ    v?-1^ 

where  T   is  the  amount  of  time  that  x(t)  falls  in  the 
range  (x,  x+dx)  during  an  observation  time  T. 

For  this  experiment,  the  sea  surface  displacement  was 
assumed  to  be  the  sum  of  an  infinite  number  of  independent 
components  of  random  phase.   Under  these  conditions  the 
measured  frequency  distribution  of  the  sea  surface  displace- 
ment tends  to  be  Gaussian.   For  zero  mean  the  Gaussian 
probability  density  function  is  given  by 

2 
PGU)  =  —  e?  /2  (5.13) 

where  £  is  equal  to  the  value  of  the  measured  quantity 
divided  by  the  standard  deviation.   However,  if  the  non- 
linearities  of  the  boundary  condition  become  important, 
departures  from  Gaussian  are  expected.   Longuet-Higgins 
(1963)  showed  for  "moderately"  non-linear  swell  in  deep 
water  that  the  probability  distribution  of  ri  is  better 
described  by  the  Gram-Charlier  probability  density  function. 
pQC(0  is  given  by: 
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2 

m^       m^      m~ 

pgc(c)  =  pGU)  [i.oi^iy  (-^  H]|  +  -^-  h6)  +  ...3.  (5.i4) 


Here  the  m  '  s  are  cumulants  with  nu  being  the  skewness  of 

the  data  record,  m^  the  kurtosis  minus  three,  and  H   the 

Hermite  polynomials  of  degree  n.   The  Hermit e  polynomials 
are  generated  from 


H 

n 


=  n   n(n-l)  £n  2  ,  n(n-l)  (n-2)  (n-3)  K  /,  l0 

c,       -j^j     2  2!  2  '   *  *  '   '     ' 


Therefore 


H3  =  53  -  3   , 


*J    /-  2 

h^  =  r  -  6r  +  3    , 


and 


H,  =  K6  -   15C4  +  ^K2   ~   15   .  (5.16) 


If  the  distribution  of  £  is  Gaussian,  the  skewness  is  zero 
and  the  kurtosis  is  equal  to  three  which  results  in  these 
and  all  higher  order  cumulants  being  zero  making 


PGC(0  =  PGU)  . 


Kinsman  (i960)  made  measurements  of  wind  generated  waves 
and  showed  their  distribution  was  only  approximately  Gaussian 
and  better  approximated  by  the  Gram-Charlier  distribution. 
He  found  the  distribution  exhibited  a  positive  skewness  and 
attributed  this  skewness  to  the  physical  fact  that  waves  have 
a  distinct  tendency  to  form  short  sharp  crests  and  long 
shallow  troughs. 
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F.   GOODNESS-OF-FIT  TEST 

The  chi-square-  goodness-of-fit  test  is  used  to  compare 
a  measured  probability  density  function  (PDF)  to  a  theoret- 
ical one.   A  statistic  with  an  approximate  chi-squared 
distribution  is  calculated  as  a  measure  of  the  discrepancy 
between  an  observed  PDF  and  the  theoretical  PDF.   The 
chi-squared  parameter  is: 

X2  =  I    1  F  X  (5.17) 

i=l      i 

where  f.  is  the  observed  frequency  of  occurence,  F.  the 
theoretical  frequency  of  occurence,  and  J  is  the  number  of 
intervals  into  which  the  frequency  distribution  is  divided. 
If  a  calculated  distribution  matches  exactly  a  theoretical 
distribution  the  chi-squared  parameter  is  zero. 

Two  methods  can  be  used  to  calculate  the  chi-squared 
parameter  for  the  goodness-of-fit  test.   In  the  first  method, 
the  £-axis  of  the  PDF  is  divided  into  J  intervals  of  equal 
size.   The  f.'s  in  this  case  represent  the  number  of  data 
points  in  each  interval  between  divisions  and  f .  varies 
across  the  distribution.   In  the  second  method,  the  PDF  is 
divided  into  J  intervals  of  equal  probability  density,  each 
containing  approximately  an  equal  number  of  data  points  in  f .  . 
Williams  (1950)  showed  that  to  give  the  test  sufficient 
power  the  number  of  intervals,  J,  should  be  given  by 


J  =  2 


\J 


2(n-^  (5.18) 


i\o 


where  n  is  the  number  of  time  samples.   C,  under  the  null 
hypothesis,  is  the  cutoff  value  for  the  level  of  significance 
which  as  used  in  this  analysis  was  taken  as  0.05. 
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VI.   ANALYSIS  OF  DATA 

The  analog  data  were  digitized  at  0.2  second  intervals. 
The  relationship  between  the  upper  limit  of  the  bandwidth 
to  be  investigated,  f  ,  and  the  sampling  interval  At,  is 


ife  ■  2lt  •  (6'X) 


Here  f   is  the  Nyquist  frequency  which  in  this  experiment 
was  established  at  2.5  Hz.   This  was  sufficiently  high  to 
avoid  aliasing  or  folding  problems.   The  maximum  lag  time 
was  chosen  as  consistent  at  10  percent  of  the  record  length. 
This  results  in  twenty  degrees  of  freedom  for  each  computed 
spectra.   The  80  percent  confidence  limits  for  twenty 
degrees  of  freedom  are  between  0.71  and  1.60  of  the  measured 
spectral  estimates  according  to  the  chi-squared  distribution 
Because  the  record  lengths  were  different  for  each  run, 
there  was  a  different  frequency  resolution  for  each  run. 
The  frequency  resolution  for  each  run  is  shown  on  the 
respective  figures. 

A.   VERTICAL  VELOCITY  MEASUREMENTS 

The  predicted  and  measured  vertical  velocity  energy- 
density  spectra,  and  the  phase  angle,  and  coherence-spectra 
of  the  wave  height  and  vertical  velocity  were  computed  for 
runs  two  through  six  at  the  depths  shown  in  Table  I. 
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1 .   Predicted  and  Measured  Vertical  Velocity 
Energy-Density  Spectra 

The  vertical  velocities  were  predicted  using  (4.4). 

The  predicted  vertical  velocity-spectrum  was  determined  from: 

|w(f)|2  =    c?i  smh  kCh+z)^!      ^2  (6    } 

11       f    smh  kh      '  ' v   ' 

The  term  in  brackets  is  the  square  of  the  transfer  function, 
(4.6),  relating  the  vertical  velocity  to  n. 

Measured  vertical  velocity  energy-density  spectra, 
for  Runs  Two  and  Six,  are  shown  in  Figures  (11)  and  (12). 
The  principle  energy  band  was  from  0.04  to  0.37  Hz  as  deter- 
mined from  the  wave  spectra.   For  Run  Six,  depth  5.70  meters, 
the  theoretical  values  of  the  vertical  velocity  spectrum 
were  within  three  to  five  percent  of  the  measured  vertical 
velocity  spectrum  values  over  a  frequency  range  of  0.04 
to  0.34  Hz.   The  true  error  of  the  velocities  is  the  square- 
root  of  the  spectral  differences  and  was  approximately  two 
percent.   This  is  within  the  combined  error  of  the  wave  and 
velocity  meters;  hence  most  of  the  measured  energy-density 
within  this  bandwidth  was  wave  induced  with  relatively 
little  spectral  energy-density  being  due  to  turbulence. 

On  the  other  hand,  for  Run  Number  Two  (the  run  of 
deepest  depth  at  14.62  meters)  the  absolute  difference 
between  the  measured  and  theoretical  spectra  was  small  over 
the  much  narrower  bandwidth  of  0.05  to  0.19  Hz.   The  large 
absolute  difference  between  the  measured  and  theoretical 
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Figure  11.   Measured  and  Theoretical  Vertical  Velocity 
Spectra  for  Run  Number  Two.   Depth  1*1.62  Meters. 
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spectra  at  high  frequencies  for  both  runs  indicates  that 
beyond  0.19  Hz  for  Run  Number  Two  and  0.3^  Hz  for  Run 
Number  Six,  that  most  of  the  energy  was  due  to  turbulence. 

The  decrease  of  wave  induced  motion  with  increased 
depth  can  also  be  seen  by  comparing  the  vertical  velocity 
spectra  for  the  two  runs.   For  Run  Number  Six  there  was  a 
peak  in  energy  due  to  a  three  second  wave  component.   This 
peak  was  not  present  in  the  spectra  for  Run  Number  Two. 
That  the  energy  present  at  this  frequency  and  depth  was 
due  to  turbulence  can  be  seen  by  the  large  deviation 
between  the  measured  and  theoretical  spectra. 

Like  the  spectral  energy-density  in  the  wave  induced 
part  of  the  spectrum,  the  turbulent  energy-density  level 
was  one  order  of  magnitude  less  for  Run  Number  Two  compared 
to  Run  Number  Six.   This  indicates  that  turbulent  motion 
decreased  at  the  same  rate  as  wave  induced  motion  under  the 
particular  conditions  of  moderate  swell,  low  winds,  and 
intermediate  water  depth  and  that  the  intensity  of  the 
turbulence  is  associated  with  the  rms  value  of  the 
wave-induced  motion. 

A  measure  of  the  total  energy  available  in  the 
measured  spectra  is  given  by  the  variance.   The  value  of 
the  variance,  and  therefore  the  total  kinetic  wave-induced 
energy  present,  decreased  with  increasing  depth.   There  is 
almost  an  order  of  magnitude  difference  between  the  two 
values  indicating  a  rapid  decrease  in  energy  with  depth 
as  predicted  by  theory. 
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2.  Coherence  of  the  Wave  Height  and  Measured 
Vertical  Velocity 

The  coherence  spectra  are  also  shown  in  Figures 
(11)  and  (12).   In  general,  the  coherence  was  0.8  or 
greater  throughout  the  range  of  significant  wave  induced 
energy-density.   This  means  that  there  is  a  high  correla- 
tion between  the  waves  and  wave-induced  vertical  motion  and 
that  relatively  little  of  the  energy  in  the  significant 
energy-density  band  was  due  to  turbulence,  but  instead, 
was  wave-induced . 

3.  Phase  Angle  Between  Wave  Height  and  the 
Vertical  Velocity 

The  phase  angle  theoretically  leads  wave  height  by 
90  degrees.   The  two  computed  phase  spectra  in  Figures  (11) 
and  (12)  further  reflect  the  accuracy  of  linear  theory  at 
both  shallow  and  deep  depths. 

Bordy  (1972)  encountered  the  greatest  departure  from 
linear  wave  theory  in  the  computation  of  the  phase  angle. 
Much  effort  was  expended  during  the  course  of  this  experi- 
ment to  ensure  that  phase  angle  errors  due  to  experimental 
technique  were  minimized.   The  concurrance  of  the  constant 
phase  angle  with  depth  further  supports  linear  theory.   It 
is  felt  that  keeping  the  wave  sensor  directly  over  the  flow- 
meter was  the  sole  reason  for  the  difference  in  the  results 
obtained  in  this  research  as  compared  with  those  of  Bordy. 

l\ .   Comparison  of  Phase  Angle  and  Coherence  for 
Different  Depths" 

As  can  be  seen  from  Figure  (13)  >  there  is  a  definite 

relation  between  the  phase  angle  for  a  given  frequency  and 
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Figure  13.   Vertical  Phase  Spectra  for  Runs  Two  Through  Six 
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the  depth  of  measurement.   The  90  degree  phase  shift  was 
observed  to  a  period  of  0.1!  Hz  at  the  shallowest  depth, 
but  only  as  far  as  0.25  Hz  for  the  deepest  case. 

The  upper  frequency  limit  of  high  coherence  also 
decreases  with  increasing  depth  as  illustrated  in  Figure 
(1*0.   Both  of  these  phenomena  can  be  explained  by  linear 
theory.   Linear  theory  predicts  decreasing  particle  motion 
with  increasing  depth  due  to  a  passing  wave  and  that  for  a 
given  depth  the  lower  frequency  waves  induce  greater  parti- 
cle motion.   This  is  reflected  in  the  coherence  and  phase 
plots.   The  highest  frequency  at  which  the  coherence 
approached  one  and  at  which  the  phase  compared  well  with 
theory  for  all  analyzed  runs  corresponded  with  the  onset 
of  significant  turbulent  energy  contribution. 

B.   HORIZONTAL  VELOCITY  MEASUREMENTS 

The  predicted  and  measured  total  horizontal  velocity 
energy-density  spectra  were  computed  for  Run  Number  One. 
The  phase  angle  and  coherence  spectra  were  computed  using 
the  wave  height  and  v-component  of  velocity. 

1.   Predicted  and  Measured  Horizontal  Velocity  Spectra 
The  total  theoretical  horizontal  velocity  was 
computed  from   (4.3).   The  theoretical  energy-density 
spectrum  of  the  total  wave-induced  horizontal  velocity  was 
calculated  from  the  expression 


■)"/„v|2    r  f   cosh  k(h+z)-,2  i  ,„,i2         ,c    ~s 
u<f)l   ■  [2¥   sinh  kh ]   ln(f)l   •       C6-3) 
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Figure  14.   Vertical  Coherence  Spectra  for  Runs  Two 
Through  Six. 
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The  term  in  brackets  is  0.5),  the  transfer  function 
relating  the  total  horizontal  velocity  to  n .   The  total 
spectrum  was  found  by  calculating  the  individual  u-  and 
v-component  spectra  and  summing  over  the  frequency  band 
of  interest.   The  spectra  are  illustrated  in  serni-log  format 
in  Figure  (15). 

The  difference  between  the  spectra  at  the  peak  was 
found  to  be  about  12  percent,  giving  an  absolute  value  of 
the  difference  of  about  3.5  percent  in  the  significant 
energy-density  frequency  band. 

2 .  Coherence  of  Wave  Heights  and  Measured  Total 
Horizontal  Velocity 

The  coherence  (of  n.  and  the  v-component  of  velocity) 
was  calculated  to  be  greater  than  0.9,  Figure  (15)  5  through- 
out the  significant  energy  frequency  band  of  0.03  to  0.07 
Hz.   This  high  coherence  was  due  to  the  fact  that  the  v- 
component  of  velocity  was  primarily  wave-induced  and  the 
turbulent  contribution  to  the  spectrum  was  small. 

3.  Phase  Angle  Between  Wave  Height  and  the  Total 
Horizontal  Velocity 

The  theoretical  phase  angle  for  the  wave  height 

and  the  v-component  of  velocity  is  zero  degrees.   This  was 

found  to  be  true  in  this  case  as  illustrated  again  in 

Figure  (15). 

C.   PROBABILITY  DENSITY  FUNCTIONS  AND  STATISTICS  OF  THE 
DATA  SETS. 

A  probability  density  function  was  computed  for  each 

measured  quantity  for  the  six  runs.   The  points  comprising 
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Figure  15.   Measured  and  Theoretical  Horizontal  Velocity 
Spectra  for  Run  Number  One.   Depth  7.04  Meters. 
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a  histogram  of  the  data  sets  were  plotted  against  both  a 
Gaussian  and  a  Gram-Charlier  distribution.   The  chi-squared 
goodness-of-f it  test  was  computed  for  both  theoretical 
distributions.   The  statistics,  variance,  standard  deviation, 
skewness,  and  kurtosis  were  also  computed  for  all  eighteen 
data  sets.   These  results  are  listed  in  Table  IV. 
1.   Probability  Density  Functions 

A  typical  PDF  is  shown  in  Figure  (16).   It  is 
known  that  the  ocean's  surface  has  approximately  a  Gaussian 
distribution.   However,  the  non-linearities  of  the  surface 
introduce  values  of  the  skewness  and  kurtosis  that  deviate 
from  Gaussian  and  result  in  a  Gram-Charlier  distribution. 
This  is  explained  by  the  physical  fact  that  waves  tend  to 
form  short,  sharp  crests  and  long  troughs.   For  the  waves 
this  produces  a  positive  skewness,  indicating  more  values 
in  the  time  series  below  mean  water  level  than  above. 

It  follows  that  the  wave  induced  particle  motion 
should  follow  the  same  behavior.   However,  of  the  eighteen 
data  sets  analyzed,  three  did  not  have  positive  skewness. 
These  three  were  the  velocity  components  u  and  w  from  Record 
Five  and  w  from  Record  Three.   Why  this  is  so,  despite  the 
fact  that  all  the  wave  records  were  positively  skewed, 
remains  unexplained.   A  possible  explanation  could  be  that 
this  peculiarity  is  caused  by  internal  waves  which  were 
prevalent  during  the  experiment. 
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2.   Goodness-of-Fit  Test 

The  chi-squared  goodness-of-f it  test  was  used  to 
compare  the  computed  PDF  to  both  the  Gaussian  and  Gram- 
Charlier  distributions  for  all  eighteen  data  sets. 
Qualitatively,  based  on  the  smaller  chi-squared  values 
obtained  from  a  comparison  of  the  PDF  with  the  Gram- 
Charlier  distribution,  it  can  be  stated  that  all  eighteen 
records  more  closely  fit  the  Gram-Charlier  distribution 
than  the  Gaussian.   The  chi-squared  parameter  for  all 
eighteen  data  sets  is  also  listed  in  Table  IV. 


56 


VII.   CONCLUSIONS 

The  measurements  in  this  research  were  made  under  condi- 
tions of  moderate  swell  and  low  wind  velocity.   Wave  height 
was  measured  with  a  penetrating  wave  staff  and  water 
particle  motion  measured  with  an  electromagnetic  flowmeter. 

A.  ENERGY-DENSITY  SPECTRA 

The  wave-induced  vertical  velocity  energy-density 
spectra  were  computed  for  six  depths.   The  frequency  of  the 
peak  energy-density  remained  constant  with  depth.   However, 
the  significant  energy-density  of  the  wave-induced  motion  at 
high  frequency  decreased  with  increasing  depth  according  to- 
the  attenuation  aspects  of  linear  theory.   The  magnitude  of 
the  energy-density  spectra  also  decreased  with  Increasing 
depth  according  to  theory.   It  was  noted  that  the  lowest 
frequency  at  which  turbulence  contributed  significantly  to 
the  spectra  shifted  to  lower  frequencies  with  increased 
depth.   The  intensity  of  the  turbulent  contribution  decreased 
with  depth  and  appears  associated  with  the  variance  of  the 
wave-induced  velocities. 

B.  COHERENCE 

Coherence  of  the  wave  height  and  wave-induced  velocities 
was  computed  to  be  greater  than  0.85  throughout  the  range 
of  significant  energy-density.   The  upper  frequency  limit  of 
high  coherence  decreased  with  increasing  depth.   The  value 
of  the  upper  limit  was  established  at  the  frequency  where 
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turbulence  became  an  important  factor  in  the  energy-density. 
The  near  unity  values  of  coherence  over  the  band  of  signifi- 
cant spectral-density  imply  a  near  linear  process  between 
wave  height  and  induced  water  particle  motion.   It  is  felt 
that  the  small  irregularities  in  coherence  in  the  significant 
energy-density  band  were  due  to  turbulence. 

C.  PHASE  ANGLE 

The  computed  phase  angle  spectra  of  the  wave  height  and 
induced  water  particle  motions  were  consistent  with  the 
theory  at  all  depths  for  both  vertical  and  total  horizontal 
velocity  cases.   This  is  the  part  of  the  analysis  where 
Bordy  (1972)  encountered  the  greatest  deviation  from  theory. 
It  is  imperative  that  the  wave  staff  be  directly  over  the 
flowmeter  in  order  that  there  be  no  phase  error.   A  similar 
decreasing  upper  frequency  limit  with  increasing  depth  for 
consistent  phase  angle  values  was  present  as  it  was  for  both 
the  coherence  and  the  energy-density  spectra. 

D.  STATISTICAL  INTERFERENCES 

The  probability  density  functions  of  all  measured 
quantities  were  qualitatively  found  to  be  better  described 
by  the  Gram-Charlier  distribution  than  by  the  Gaussian  using 
the  chi-squared  goodness-of-fit  test.   The  variance  of  the 
vertical  velocities,  a  measure  of  the  kinetic  energy  present, 
decreased  with  increasing  depth  which  is  consistent  with 
linear  theory.   The  skewness  of  the  wave  height  frequency 
distribution  was  found  to  be  positive  in  all  cases  as  a 
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consequence  of  the  fact  that  waves  tend  to  have  short, 
sharp  crests  and  elongated  troughs.   The  values  of  skewness, 
with  few  exceptions,  were  found  to  be  positive  for  the  wave- 
induced  quantities  as  well.   These  exceptions  might  be 
explained  by  the  presence  of  internal  waves. 
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